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We calculate the one-particle density of states for the Mott-Huhhard insu- 
lating phase of the Hubbard model on a Bethe lattice in the limit of infinite 
coordination number. We employ the Kato-Takahashi perturbation theory 
around the strong- coupling limit to derive the Green function. We show that 
the Green function for the lower Hubbard band can be expressed in terms 
of polynomials in the bare hole-hopping operator. We check our technique 
against the exact solution of the Falicov-Kimball model and give explicit re- 
sults up to and including second order in the inverse Hubbard interaction. 
Our results provide a stringent test for analytical and numerical investiga- 
tions of the Mott-Hubbard insulator and the Mott-Hubbard transition within 
the dynamical mean-field theory. We find that the Hubbard-UI approxima- 
tion is not satisfactory beyond lowest order, but the local-moment approach 
provides a very good description of the Mott-Hubbard insulator at strong 
coupling. 

PACS numbers: ll.lOFd, 11.21. +a, 11.30+h 



1. INTRODUCTION 

As emphasized by Mott early on0'i, interacting electrons in a single half- 
filled band may undergo a quantum phase transition from a metallic to an 
insulating state without breaking the translational or spin symmetry of the 
underlying Hamiltonian (Mott-Hubbard transition!). The explicit formula- 
tion of the corresponding minimal Hamiltonian by Hubbardcl, GutzwilleiH, 
and Kanamorii, allowed a concise discussion of the Mott-Hubbard transi- 
tion. In the Hubbard model, the electrons' on-site interaction of strength U 
competes with their kinetic energy, characterized by the bandwidth W for 
vanishing interactions. 



* Dedicated to Peter Wolfle on the occasion of his 60th birthday. 
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The Hubbard Hamiltonian poses a most difficult many-body problem, 
and approximate treatments had to be devised. For example, within the 
Gutzwiller variational schemeO Brinkman and RiceE corroborated Mott's 
view of the Mott-Hubbard transition as a continuous quantum phase transi- 
tion: At some finite critical interaction strength, Uc = 0{W), the quasipar- 
ticle weight goes to zero comine from the metallic side, and the gap closes 
coming from the insulating sideB. Qualitatively the same picture is found 
from the exact solution of the 1/r Hubbard model in one dimensionEru. 

A few years ago, these views were challenged by Georges and Kotliar 
and their collaborators; for a review, see Ref. They developed further 
and studied the dynamical mean-field theory for the Hubbard model which 
becomes exact in the limit of infinite lattice coordination numberS. Results 
from various analytical and numerical techniques led them to the conclusion 
that the Mott-Hubbard transition is discontinuous in the sense that the 
gap jumps to a finite value when the quasiparticle weight becomes zero at 
f/c,2 = Uc (scenario of a 'preformed gap'). They argue that the insulating 
solution to the equations persists down to Uc,i < C/c,2 but its energy is higher 
than that of the metallic state. 

Various approaches to the Mott-Hubbard transition for lattices with 
infinite coordination number yield results in favor of either of the conflicting 
scenarios; for a review, see Refs. ^J^, and Refs. 11-17 for more recent treat- 
ments. Insightful physical arguments helped to sharpen the physical and 
mathematical implications of the scenario of a preformed gap transition^^"^^ 
but cannot resolve the issue. Since an exact solution to the dynamical mean- 
field equations is still lacking, all analytical approaches^^"^^'^"^ are neces- 
sarily approximate in nature. Numerical investigations of the dynamical 
mean-field equations^^"^'''^^''^^ involve (i) discretization, (ii) numerical diag- 
onalization, (iii) interpolation, (iv) iteration of the self-consistency cycle, 
and (v) extrapolation to the thermodynamic limit. The Random-Dispersion 
Approximation0 provides an alternative to the dynamical mean-field ap- 
proach. However, apart from (iv), similar steps are required in its numerical 
implementation. Therefore, the region of applicability of the available nu- 
merical techniques is not a priori clear either. 

In such a situation, benchmark tests need to be provided in order to 
assess the quality of the various analytical and numerical approaches deep 
in the metallic and insulating phases where perturbation theory provides 
reliable answers. If an approximate technique fails in the perturbatively ac- 
cessible regimes where others may pass, its predictions on the Mott-Hubbard 
transition should not be taken too seriously. 

In this work, we provide such a benchmark test for the Mott-Hubbard 
insulator at zero temperature. In Sect. |2] we introduce the definitions of 



Mott Hubbard Insulator in Infinite Dimensions 



the Hubbard model and the Fahcov-Kimball modelEU and the relevant one- 
particle Green functions. We restrict ourselves to the half-filled Bethe lat- 
tice with an infinite coordination number, and exclude all possible ordered 
ground states which would conceal the Mott-Hubbard transition. 

The ground state of the Mott-Hubbard insulator is vastly degenerate. 
Therefore, in Sect. we invoke Kato-Takahashi perturbation theory to 
calculate the one-particle Green function systematically in the insulating 
regime. As a simple application, we re-derive the well-known lowest-order 
solution^!. 

In Sect. we use the exactly solvable Falicov-Kimball model as a test 
caseEH. We show how to overcome the resummation problem of our series 
expansion. We find excellent agreement between our results to second order 
in \/U and the exact density of states down to C/ ~ 0.75VF; for the Falicov- 
Kimball model, = W/2 for the transition. 

This observation makes us confident that our second-order results for 
the Hubbard model. Sect. [5], provide a reliable test down toU ^ 1.5 W . Our 
approach allows an estimate of the critical value for the closure of the gap. 
This value corresponds to Ud in the scenario of a discontinuous transition, 
and to Uc if the scenario of a continuous quantum phase transition is the 
correct one; for the Hubbard model, there is consensus that the insulator no 
longer exists below 0.91^ < Uc,i < 1.3W. 

In Sect. ^ we make a first comparison with two other analytical meth- 
ods. The Hubbard-III approximation turns out to be inappropriate whereas 
the local-moment approach provides a very good description of the Mott- 
Hubbard insulator. A more detailed comparison is planned to be published 
elsewhere. Short conclusions. Sect. |7], close our presentation. 



2. DEFINITIONS 
2.1. Hamilton Operators 



In this work we study the Hubbard mo^ 




H = f + Ub . (1) 
The kinetic energy of the electrons reads 

^ = ^ E <.^^+r,. , (2) 

where cf^, ci^a are creation and annihilation operators for electrons with 
spin a at site z, where i runs over all lattice sites and r denotes the Z nearest- 
neighbor vectors. Since we are interested in the Mott insulating phase, we 
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consider exclusively a half-filled band where the number of electrons N equals 
the number of lattice sites L. 

We restrict our analysis to the Bethe lattice where each site is connected 
to Z other sites without the generation of loops, and the limit Z — > oo is 
implicitly understood henceforth. Then, the non-interacting density of states 
becomesEy 

where W = At is the bandwidth. In the following, we shall set i = 1 as our 
energy unit. 

The Hubbard interaction reads 

°=i:("«4)("-'4) • « 

where nj^o- = cf^Ci^a are the density operators for a given spin a. For later 
use we also define the operators for the local electron density, fii = hi^-\ +hi^^^. 
The operators Sf = (c+^^Q^j^ + cf^^Ci^f)/2, Sf = (c+^q^i - c+j^Ci,|)/(2i), and 
Sf = (rij^i — nj^j)/2 are the three components of the spin-1/2 vector operator 

We have chosen the chemical potential in such a way that the Hamil- 
tonian explicitly exhibits a particle-hole symmetry, i.e., = guarantees 
a half- filled band. We restrict our analysis to the Mott-Hubbard insulator 
at zero temperature which is characterized by a finite ground-state entropy 
densityi'i, s = ln(2) + C'(ln(L)/L). This means that in the Mott-Hubbard 
insulator each lattice site is equally likely occupied by an electron with spin 
I or I, irrespective of the spin at any other lattice site. 

Later we shall also address the case of the Falicov-Kimball model0 (or 
simplified Hubbard model) in which only one species of electrons is mobile, 

^"^ = 4e^m+.. (5) 

The mobile c electrons interact with the immobile / electrons via the Hub- 
bard interaction, 

= E kc. - \) (ftn - 1) . (6) 



The immobile electrons are randomly distributed over the lattice such that 
each lattice site is occupied by an / electron with probability p = 1/2, 
irrespective of the / occupation at any other lattice site. 
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2.2. Green Functions 

The time-dependent local single-particle Green function at zero temper- 
ature is given byEiJ 

G{t) = -■^\T.^mAt)ct]) ■ (7) 

Here, T is the time-ordering operator, (...) implies the average over all 
ground states with energy £"0, and {fi = 1) 

Ci,a{t) = exp{iHt)ci^a exp{-iHt) (8) 

is the annihilation operator in the Heisenberg picture. Recall that the sym- 
metries of H are unbroken so that each lattice site gives the same contribu- 
tion in (|^. 

The Bethe lattice is bipartite so that we may define (—1)* = +1 for the 
A sites of the lattice which are surrounded by B sites only (and vice versa) , 
for which (—1)* = — 1 {i € B). Then, the particle-hole transformation 

'T: c+.-(-irQ,. ; Q,.^(-irc+, (9) 

maps H onto itself. We can readily identify the contributions from the lower 
(LHB) and upper (UHB) Hubbard bands to the Fourier transform of the 
local Green function (77 = 0"^), 

/oo 
dte'^'G{t) = Glub{oj) + Guhb(c^) , 
-00 

Glhb(c^) = jJ2(^ta\^ + iH-Eo)-ir]\~'ci,^) , (10) 



G'uhb('^) — — Glhb(— w) 5 

due to particle-hole symmetry. Therefore, it is sufficient to evaluate the local 
Green function for the lower Hubbard band which describes the dynamics 
of a hole inserted into the ground state. 

The density of states for the lower Hubbard band can be obtainejtLfrom 
the imaginary part of the Green function ([lO| ) for real arguments viaE^ 

-Dlhb('^) = -9{Glhb('^)} , /^Ehb < ^ < /^LHB < • (H) 
vr 

Due to particle-hole symmetry we have D\jiib{uj) = Dlhb(— i^)- The gap in 
the Mott-Hubbard is symmetric around a; = so that the one-particle gap 
in the Mott-Hubbard insulator is given by 

A = 2|/i+HBl >0. (12) 
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For the Falicov-Kimball model we shall address the local Green function of 
the mobile electrons for the lower Hubbard band 

Gl^Bi^) = ^ E {^t [u^ + iH- Eo) - iv] £i) • (13) 

Again, the local Green function is given by G^^(aj) = Gi^^^{lo) — Gf^^{—uj) , 
the density of states is symmetric around to = 0, and Eq. ( [Ill ) hold analo- 
gously. For later reference we further define 



3. STRONG-COUPLING EXPANSION 

3.1. Kato Takahashi Perturbation Theory 

Based on Kato's degenerate perturbation theoryU, TakahashiH devel- 
oped the perturbation expansion in 1/U for the Hubbard model at zero 
temperature. For large interaction strengths U, the ground states \ipn) of H 
in (|^) can be obtained from states \(j)n) without a double occupancy, 

IV'n) = f |<^n) , Po\<Pn) = \<Pn) , (15) 

where Pj projects onto the subspace with j double occupancies. 

Eq. ( [Tsl ) is readily interpreted. In the large-coupling limit, T in (|l]) 
is considered as a perturbation on UD. Addressing the ground states we 
therefore start from eigenstates with zero double occupancies into which 
the operator F successively introduces at most m double occupancies to 
generate the ground states IV'n) of H to m-th order in t/U. The operator 
r reduces all operators to the subspace with zero double ocpipancies. In 
particular, F+F = Pq so that overlap matrix elements obeyH {ijjmli^n) = 

{(pm\<t>n)- 

The Schrodinger equation for the ground states 

Ht\(Pn) = ^of (16) 

leads to 

h\(Pn)=0 , h = t+Ht-EoPo , (17) 

i.e., the eigenvalue does not change under the transformation. The creation 
and annihilation operators are transformed accordingly, 

^ = r"^C^^f , ^ Cl^^ = f +Q,„f . (18) 
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The derivation of the exphcit expression for F can be found in Refs. 32,33, 
where it is shown that 

-1/2 



Here, 



f = PPo (PoM) ■ (19) 



P = Po-J2 J2 ^'''^ • • • TS''"+' , 

n=l ri + ...+rn+i=n 

5° = -P\ (20) 

TTr Z-^ Ar 

The term to n-th order in P contains all possible electron transfers generated 
by n applications of the perturbation T; its contribution is proportional to 
(i/C/)". The square-root factor in (|l9|) guarantees the size-consistency of the 
expansion, i.e., it eliminates the 'unconnected' diagrams in a diagrammatic 
formulation of the theoryll. The square root of an operator is understood 
in terms of its series expansion, i.e.. 



. . . ^(2n-l)!!r- - - in 

(Wo) ^^o + E ^ (2^),, [Po{Po-P)Po] . (21) 



The Green function for the lower Hubbard band becomes 



Glhb(w) = -E(4^ f^ + ^-i^l ^ Ci^a) , (22) 



L 

where (...) now implies the average over all states \(j)n) with no holes and no 
double occupancies (A^ = L electrons), i.e., the average over all single-spin 
configurations. 

Up to and including second order in 1/U a straightforward expansion 
gives 

f = Po + S^fpQ - S^fPofPo + S^fS^fPo - ^Pofs'^fPo , (23) 

Q,. = PoC^,aPo + PoCi,J^fPo + P^f f Pq + P^h^J^fS^fP^ 

-^PoCi,^PofS^fPo - ]^P^fS^fP^c,,„Po . (24) 

In the derivation of (^) we used the fact that Cj^o- in ( p2D acts on states \(j)n) 
with no holes and no double occupancies. The transformed Hamiltonian up 
to second order reads 

h = ^{L-2N) + hQ + hi + h2 + ... (25) 
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with 



ho 
hi 

h2 



PoTPo 



L 



PofS^fPo + — , 
PofS^fS^fPo - ^PofPofS^fPo - ^PofS^fPofPo 



(26) 



where we used E^^ = 0, E^^ = —L/(2U), and E^^^ = for spin-disordered 
ground stated. 

Including all second-order terms, the Green function for the lower Hub- 
bard band ( p!o[) becomes 



G2{^) = jJ2 \4,a \^ + + hQ + hi + h2-\r] 



-1 



(27) 



For the Falicov-Kimball model, all operators T need to be replaced by 
from (^), and the spin index in ( p4|) and ( p7|) must be dropped. 



3.2. Green Function to Lowest Order 

As a first application we derive the well-known lowest-order result for 
the local hole Green function, 

G^{z) = i ^ /c+, [z + ho] ' , (28) 

where 2: is a complex variable which we shall put to z = uj + U/2 — i-q at the 
end of the calculation. 

A series expansion for Go{z) in terms of ho/z converges for \z\ > 2, 

The term in n-th order describes the creation of a hole at site i which moves 
through the lattice in such a way that it returns after n nearest-neighbor 
hops. Therefore, n is even, and the spin background is restored after the 
hole has returned to i. 

Obviously, the hole can return to i more than once during its excursions. 
We define S{z) as the contributions to the series such that the hole does not 
return in intermediate steps. Then, 

1 °° 1 
G„W = -i:[5(.)r = j|^-^. (30) 
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On a Bethe lattice there are no loops. Therefore, the excursions which 
contribute to S{z) start with a jump of the hole to one of the Z nearest 
neighbors i + t which serves as a starting point for arbitrary excursions 
which return to i + r and avoid i. The latter constraint is irrelevant for 
Z — > oo. The last jump takes the hole from site i + r back to i. This implies 

2 oo 



S{z) = Z 



1 



z^/Z 



m=0 



1 



z^{l-S{z)) ■ 



This quadratic equation is readily solved to give 

1 

2 



Go{z) = zSiz) 



(31) 



(32) 



The sign has to be chosen to fulfill 



in agreement with the exact solutior 
the proper limit for ?R.{z} — oo. 

We perform an analytical continuation from \z\ > 2 to the whole com- 
plex plane so that the formula (|3^ ) is valid for all z. In particular, for 
z = uj + U/2 — IT] the density of states of the lower Hubbard band becomes 

1 

2^ 

(33) 



/4 - (w + [7/2)2 
p{uj + U/2) . 



U/2\ < 2 



This result does not come as a surprise because the motion of the hole 
restores the spin configuration. Therefore, the local density of states of the 
non-interacting system (P) is reproduced apart from the energy shift U/2. 
The one-particle gap A ([l^) becomes 



Ao(C/) = [/-4 



(34) 



to leading order in l/U because /Ulhb — ~^ /2 + 2. 

For the Falicov-Kimball model we proceed along the same lines. Note, 
however, that (i) the creation of a hole at site i has only probability p = 1/2, 
cf. (pO|), and, (ii), the probability for a hole hop from i to i-|-r is also reduced 



(35) 



by a factor of p, cf. (|3l|). Thus, 



1 1 



n=0 



2z(l - SFK(z)) 



and 



5^K(z) 
Therefore, 



-Z 



zVZ 



2 oo 



m=0 



2z^l-S^^{z)) 



zS'^'^iz) 



+ ^Jz^ 



(36) 



(37) 
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in agreement with the exact solutionc3; see Sect. 4.1 



After the analytical continuation to all z and the replacement z = w + 
U /2 — \'q, we find for the density of states of the lower Hubbard band in the 
Falicov-Kimball model 



/^(a; + C//2) 
1 

2^ 



'2- (c^ + [//2)2 
The one-particle gap A (12) becomes 

to leading order in \/U because [/iLHB]"'' ~ 



+ U/2\ < \/2 . (38) 



2V2 

-U/2 + V2. 



(39) 



4. FALICOV KIMBALL MODEL 

The calculation of the higher-order terms of the Green function leads to 
a resummation problem. As its solution we propose the scheme of an effective 
operator for the hole motion. The validity of our approach is checked against 
the exact solution of the Falicov-Kimball modeled. We shall drop the index 
'FK' in this section. 



4.1. Exact One-Particle Green Function 



As shown in Ref. |29|, the one-particle Green function for all U is obtained 
from the solution of the cubic equatiorS 

[G{uj)f - 2uj[G{uj)f + G{uj){l +u}^ - U^A) - cu = . (40) 

The gap is explicitly given byil 

1/2 



A{U) = 2 



1 5^72 C/4 
■4 + — + T 



' 1 3C/2 [/6 
16 ^ ~ ^ ~ ^ ~2' 



(41) 



for U > = 2. The 1/U expansion of the gap to first and second order 
reads 



Ai(C/) = ?7- 2^/2 + 1, 



(42) 
(43) 
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All higher-order corrections turn out to be positive. It is therefore seen that 
the Hubbard bands tend to 'repel' each other in the Falicov-Kimball model 
so that the gap bends upwards when the interaction is reduced towards the 
transition. 



4.0 



I ' I ' — 

Exact result 

1/U-expansion 



3.0 - 



5. 2.0- 
<l 



1.0 - 



0.0. 



Fig. 1. Mott-Hubbard gap for one-particle excitations in the Falicov- 
Kimball model. The exact result (dashed line) and the result from the 
expansion to second order in \/U (solid line) are shown. 



The large-C/ expansion can also be used to make predictions about the 
value of the critical interaction strength Uc for the Mott-Hubbard transition. 
From A„(C/) = (n = 0, 1, 2) in (|3|), (H|), (||) we find 



C/i°) = 2.83 , up = 2.41 , up = 2.24 , ... , C/c = 2 . (44) 

It is seen that the second-order result is only about 10 percent off the exact 
result. In Fig. [l| we compare the result to second order ( ^3|) with the exact 
result (pT[). The agreement is very good for f7 > 3 = (3/2){7c- As we shall 



show in Sect. 4.3.3. the \/U expansion works equally well for the density of 
states. 
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We now set u; = z — C//2 in ( |40| ) and consider \z\ = 0(1) for the lower 
Hubbard band. This leads to 

[G{z)f-2z[G{z)f+G{z){l+z^)-z+U [[G{z)f - zG{z) + I/2] = , (45) 

from which the solution to lowest order in 1/U immediately follows as Go (-2), 
Eq. 



4.2. Green Function to First Order 

The first-order corrections Gi (w) stem from two sources, (i) the expan- 
sion of the Fermi operators (^), which are henceforth called shape-correction 
terms, and, (ii), the expansion of the Hamiltonian (p5[). The latter terms 
will lead to the gap renormalization problem. 



4-2.1. Shape- correction Terms 



We start from (|27|). The following first-order term arises due to the 
expansion (^^, 



GiA^) = jYl (rs'ctPo [z + ho 



(46) 



together with its Hermitian conjugate. The operator Cj in (46) can only act 
if there is no immobile / electron at site i. On the other hand, the operator 
jcf requires an immobile / electron at i. Therefore, the first-order shape- 
correction terms vanish for the Falicov-Kimball model. 



4-2.2. Gap Renormalization Problem 
To first order in 1/U we are left with the calculation of 



Gi(z) 



G^(z) 



^(c+ \z + ho + hi 



1 

L 



Goiz) + Giiz) + O (U 



Ci 



1 



E 



Z^hr 



h^ 



Z + hr 



Ci 



(47) 



(48) 



with z = ijj + U /2 — ir] at the end of the calculation. 
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The operator hi in (p6[) can be decomposed into two-site and three-site 
contributions!, hi = hf + hf^, 

= ~JU^' ^ c+ ,^Q/+/,c+Q+,,Po . (49) 

The three-site terms do not contribute to Gi{z): The operator hf^ allows 
the hole at I + T2 to jump over the immobile / electron at site / to reach site 
l + Ti. However, the site / cannot be passed again by the simple hole motion 
via /iQ. Therefore, the hole, created at site i, cannot return. 

The two-site terms give rise to the constant 1/{2U) when applied to a 
state with a single hole. Thus, 

It is evident from ( [50| ) that the first-order term diverges at u; = —U/2+\/2 = 
[a*lhb]o^) where the density of states to lowest order ( |3^ ) becomes zero. This 
had to be expected because the density of states and the gap are indeed 
renormalized to first order in 1/U, see (]4^). As is well known from standard 
perturbation theory in U, a shift of the one-particle resonances cannot be 
generated by the calculation of finite-order corrections to the Green function 
but requires an appropriate summation of an infinite series as performed, 
e.g., via the Dyson equationEJ. For our case of a degenerate ground state, 
there are no standard means of carrying out such a resummation. In the 
following paragraph we present the solution for the insulating phase of the 
Falicov-Kimball model on the Z — > 00 Bethe lattice in terms of an effective 
operator for the hole motion. 



-2 



1 

4C7 



1 + 



(50) 



4-2.3. Effective Operator for the Hole Motion 

For a Mott-Hubbard insulator the motion of the hole through the sys- 
tem does not alter the positions of the immobile electrons. Moreover, the 
background does not have dynamics of its own. To lowest order, the hole 
dynamics is governed by /iq, and we expect that higher orders in the 1/U ex- 
pansion correspond to a modified hole hopping. These modifications can be 
expressed in terms of an effective operator for the hole motion. In fact, by 
comparing ( ^8|) and (^), the effective hopping to first order must take the 
form 

hi^hf = ^. (51) 
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We insert this result into ([47|) and find 

G'i(z) = Go(z + l/(2C/)) . 
The first-order contribution is merely a band shift by \/{2U). 



(52) 



4-2. 4- Density of States 

The density of states for the lower Hubbard band is simply shifted 
downwards in energy by 1/(2C/), 



2tt 



U 1 ' 
2 +2Z7. 





U 1 




UJ-\ \ 




2 2U 



Therefore, the gap to first order becomes 



Ai{U) = U-2V2 + ^ , 



< V2 . (53) 



(54) 



because [/Ulhb]!" = —U/2 + \/2 — 1/(2C/), up to and including first order 
in 1/U. Eq. (|5^) agrees with the exact result (^2|). 



4.3. Green Function to Second Order 

4-3.1. Shape- correction Terms 



We again start from (^7\j. The following second-order terms arise due 
to the expansion (p^, 



G2,a(z) 

G'v(^) 



i 



z + hi 



z + ho 
1 



-1 



z + h(] 



Z + hn 



-1 



PofS^c^S^f 



-^PoCiPofS^fPo - ^PofS^fPoCiP, 



YT.{TS'ctPo 



Z + hr 



-1 



(55) 
(56) 



(57) 
(58) 



together with their Hermitian conjugates. With the help of the same ar- 



guments used in Sect. |4.2j , it is not difficult to show that G2,a{z) = and 
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G2,p{z) = 0. Moreover, the series expansion as exemplified in Sect. 3.2. leads 
to G2,-y(-z) = —G2,5{z) so that the shape-correction terms vanish altogether 
for the Falicov-Kimball model also in second order in l/U. 



4-3.2. Effective Operator for the Hole Motion 



To second order in l/U we are left with the calculation of 



G2{z) 



z + hn + hi + h 



-1 



(59) 



= Go{z) + Gi{z) + G2{z) + O{u-^) 
which requires the evaluation of 



G2{z) 



L ^ 



z + ho 



{hf + hf) 



Z + hr 



z + h 



(hf + /if ) - h2 

(60) 



with z = LV + U/2 — ir] at the end of the calculation. 

With the help of the series expansion, see Sect. \i.2\ , the result for G2{z) 
can be cast into the form 



U^G2{z) 



1 1 

IT 



z + h( 



z + hr 



Go{z)-Y.{ct \z + ho 



(61) 



The first term arises from the combination of the two factors /if, the second 
term stems from the two factors /if (mixed terms are zero), and the third 
term comes from /12 in (|6C|). We note that 



L ^ 



z + hr 



c.) = -2Go(^)-5:(c^ 



z + hr 



SO that the replacements ( [5l| ) and 



h2 ^ hf 



1 



4C/2 



hn 



cA , (62) 



(63) 



in (|59D reproduce (pll). Up to and including second order we may therefore 
write 



1 



= jJ2{^t z + ho + 1/{2U) - /io/(4C/2) 



-1 



(64) 



The second-order corrections result in a shift and a narrowing of the Hubbard 
bands. 
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Fig. 2. Density of states of the lower Hubbard band in the Fahcov-Kimbah 
model, shifted by f//2, for [/ = 3 = l.bUc- The exact result (dashed line) 
and the result from the expansion to second order in 1/U (solid line) are 
shown. 



4-3.3. Density of States 

Since we know the local density of states of ho , eq. (|l4|) , the density of 
states of the lower Hubbard band becomes 



rOD f U 1 

D2{uj) = j^^dep{e)5[uj + - + €+- 



2U 



l-l/(4[/2)' 



cj + ^/2 + l/(2^) \ 
1 - l/(4[/2) ) 



(65) 



which explicitly reads 
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up to and including second order in 1/U . Therefore, the gap to second order 
becomes 

A2([/) = C/-2^+i + ^, (67) 

because 

[Mlhb]2" = + ^ - l/(2?7) - ^/2/(4?72) up to and including 

second order in \/U . Eq. ( |67D agrees with the exact result (^). 

As seen from Fig. ^, there is hardly any noticeable difference between 
the exact density of states from (p^, (40), and the second-order result (|66|). 



The favorable comparison shows that our 1/U expansion works indeed, and 
that is reliable down to Uc/U ~ 2/3 in the Falicov-Kimball model. We are 
confident that the convergence will be similarly good for the Hubbard model. 



5. HUBBARD MODEL 



In Sect. |4] we have shown that, as far as the Green function for the 
lower Hubbard band is concerned, the operator h in (|T7| ) can be replaced 
by an effective operator for the hole motion, ff^, which can be written in 
terms of a polynomial in /iq. This was possible because the / electrons are 
immobile and randomly distributed over the lattice in the Mott-Hubbard 
insulating phase so that their are no spatial correlations between the mobile 
and immobile electrons. 

Qualitatively the same applies to the Mott-Hubbard insulator in the 
Hubbard model: the spin degrees of freedom are not dynamic and spatial 
correlations between the electron spins are absent. Therefore, the replace- 
ment of the operator h by an effective operator for the hole motion, /i*^^, is 
permissible as we have verified from the explicit series expansion as used in 
Sect. 



5.1. Green Function to First Order 



5.1.1. Shape- correction Terms 



We start from (|27|). The following first-order term arises due to the 
expansion (p3), 



(68) 
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together with its Hermitian conjugate. The operator TS^ involves one ad- 
ditional hole transfer so that necessarily has the form 



GiJz) 



1 

I 



Z + hr 



-1 



(69) 



Now that we only need to determine a single real parameter si, we expand 
the right-hand-sides of (65) and (69) in 1/z and compare the terms to order 

11 , , 1 , 

(70) 



i,tj 



which reduces to 



1 

2U 



-si . 



(71) 



Therefore, the operator for the shape-correction term to first order is given 

by 



Co, 



Pn 



(72) 



5.1.2. Effective Operator for the Hole Motion 
To first order in 1/C/ we are left with the calculation of 



-1 



= Go{z) + Gi{z) + O (U-^) 
which requires the evaluation of 



Gi{z) 



L ^ 



EC 



z + hr 



z + hr 



(73) 



(74) 



with z = ijo + U /2 — irj at the end of the calculation. 

As usual, the operator hi in ( [2^ ) can be decomposed into two-site and 
three-site contributionsu, hi = hf^ + hf^, 



2s 



hi 



hf 



Lt 



SiSi+r - ^ {nm+T - 1) 



Pn 



1 

~ZU 



PO >^aa'ct^r2,aCt-A-a'Cl+ru<7'P0 , (75) 
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where Ao-o-' = 1) for a = a' {a = —a'). The operator hf^ does not move 
the hole whereas the operator hf^ moves the hole by two sites. Therefore, 
the replacement 



Z + hr 



-2 



(76) 



will lead to the desired result. The real parameters gi-2 and gi-Q are deter- 
mined as before. We expand the right-hand-sides of (|7^) and ([76|) in l/z 
and compare the terms to order 



1 1 

1 H = jY^ {eta {ai-Ahof + 5l;o) Q.a 



L 



(77) 



and to order 



(78) 



The terms are readily calculated to give 



1 

U 
3 

2C7 



5l;2 + 5l;0 , 
3(29l;2 +gi;0) 



(79) 

(80) 



where we repeatedly made use of the fact that there are no spin correlations 
in the Mott-Hubbard insulator. The solution of this set of equations gives 
= — 1/(2C/) and gi-fl = 3/{2U) so that the replacement 



hi hi 



cflF 



(81) 



is valid in (|7^). Combining (^) into (^) and (^), the Green function up 
to and including first order can be written as 



u 



z + hQ + 



3 - {h^f 
2U 



-1 



(82) 
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5.1.3. Density of States 

Since we know the local density of states of /iq , eq. (^3|) , the density of 
states of the lower Hubbard band becomes 



depie) 1 



U 



P 

1 

2^ 



U - V2U^ + 2Uuj + 3 



4- (U- \/2C/2 + 2Uio + 3 



1/2 



U 1 
2 -2Z7 



(83) 
< 2 , 



up to and including first order in 1/U. Therefore, the gap to first order 
becomes 

Ai(^) = t/-4-l, (84) 
because /^lhb ~ —U/2 + 2 + 1/{2U), up to and including first order in 1/U. 



5.2. Green Function to Second Order 

5.2.1. Shape- correction Terms 



We again start from (pT]). The following second-order terms arise due 
to the expansion (p^, 



G2,a(z) 



G'2,7(z) 



Z + hr 



-1 ^ 



Z + hr 



-1 



4 51 (^t'^ (s2;3(/io)^ + S2;l/io) + ^ ^ Ci,^\ , (85) 



G2d^) = lY.{^th\z + ho] PoCi,J^fS^f 



7 \^ta {s2;2A^of + S2;0,/3) [z + ho 



(86) 



L 



z + ho 



PofS^c,J'f 



-^PoCi,„PofS^fPo - ^PofS^fPoc^,aPo 

-1 



(87) 
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-1 » 



1 



z + h, 



Z + hr 



together with the Hermitian conjugates in (p5D-(p^). As in the previous 
section, an expansion in \/z and a comparison of the two leading orders 
fixes the unknown real parameters S2-r,i/- From (|85|)-(|88|) it follows that 



=2(2S2;3 + S2;l) 



1 1 

1 

IP 



= S2;2,/3 + S2-0,I3 
S2;2,7 + S2;0,7 



S2;2,5 + S2;0,<5 



^ = 4(5S2;3 + 2S2;1 



4f72 - 2s2;2;/3 + S2;0,/3 



(89) 

(90) 

4 f/2 - 2S2;2;7 + •S2;0,7 ' (91) 
=2S2;2;5 + S2;0,5 ■ (92) 



1 1 

1 1 



Therefore, J7^S2;3 = 1/4 and U'^S2:i = —3/4 which is also obtained if we 
use (|72D and (^) directly in ( |85|) so that the first-order shape correction term 
remains unaltered. Moreover, U'^S2-2,f3 = 1/4, C^^S2;0,/3 = —1/4; U'^S2;2,-y = 
1/4, ;72s2;0,7 = -3/4; and C/^S2;2,5 = 1/4, f/^S2;0,5 = 3/4. Including the 
contribution from the Hermitian conjugates, we find S2,2 = 2(s2;2,/3 + S2;2,7) + 
S2;2,5 = 5/(4C/2) and S2,o = 2(s2;0,/3 + S2;0,7) + •52;0,5 = -5/{iU^). Therefore, 



. 1 J 5 



so that 



Ci.( 



r-3 



(93) 



(94) 



is the shape-correction term up to and including second order in l/U. 



5.2.2. Effective Operator for the Hole Motion 
To second order in l/U we are left with the calculation of 



1 



-1 



z + ho + hi+ h2 
Go{z) + Giiz) + G2{z) + O (U-^ 



(95) 
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G2{z) 



L 





-1 






-1 . 


z + ho 




hi 


z + ho 


hi - /i2 



[z + ho 



C-i cr 



(96) 



The contribution to G2{z) which involves two factors hi can be obtained by 
the replacement hi — > (|8^) so that this term does not give a contribution 
to h2^ . For the second contribution we write 



X! i^to- z + ho h2[z + ho 



(97) 



;^E(C(^^2;3(/l0)^+52;l/l0 



Z + hr 



-2 



The equations to 0(1/ z^) and 0{l/z^) result in 



3 1 

2 1J2~ ^2;3 + 92;1 



9 1 

■4-^ = %2;3 + 252;l , 



so that 5f2;3 = 3/(4C/2), ^2;! = "S/C/^, and 

h2 ^ = ^ho {{ho)' 



(98) 



(99) 



is valid in (|9^). Combining (pS|) into (^) and (^), the Green function up 
to and including second order can be written as 



G2{z) 



1 



L ^ 



ho , 5[(/io)2-l] 



(100) 



■ . 3 - (hof ^ Sho 
z + ho-\ — h 



2U 



4C/2 



{{ho 



5.2.3. Density of States 



Up to and including second order in 1/U the density of states of the 
lower Hubbard band becomes 



D2{UJ) 



J U 3-e' 
6 luj + — + e + 



1) 



+ 



3e(e2 - 4) 



P{^u 



2U ' 4C/2 

l-ejU + 5{el-l)/{4U^) 
'l-e^/t/ + 3(3eS-4)/(4[/2) ' 



U 



2U 



(101) 
< 2 , 
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where e^o is the real sohition of the cubic equation 



U 3 -el 3e^(e2-4) , , 

Therefore, the gap to second order becomes 

A2(C/) = ^7-4-1 = Ai(?7) , (103) 

because /^lhe ~ —U/2 + 2 + l/{2U) does not change to second order in l/U. 
All even orders l/f/^" (n > 1) in the l/U expansion of the gap are zero 
for the one-dimensional Hubbard modeS, and the same might be true in 
infinite dimensions. 

Eq. ( |103D shows that the Hubbard bands tend to 'attract' each other 
in the Hubbard model. In contrast to the Falicov-Kimball model, the gap 
bends downwards when the interaction is reduced towards the transition. We 
use the large-C/ expansion to estimate the value of the critical interaction 
strength Uc for the Mott-Hubbard transition. From A„(f7) = (n = 0, 1, 2) 
in (|10|) we find 

up = 4 , up = 4.24 , up = 4.24 , ... . (104) 

The shift from uP = W to uP = 1.06W to second order in l/U is only 
six percent. We are presently carrying out the expansion for the gap to 
third order in l/U but we do not expect sizable corrections, i.e., Uc = 
(1.05 lb 0.05) is our prediction of the critical interaction strength for the 
closure of the gap. 



6. COMPARISON WITH OTHER ANALYTICAL 

APPROACHES 

6.1. Hubbard-III Approximation 

Hubbardi§ provided an early approximation for the Mott-Hubbard in- 
sulator. Its simplest form, the 'alloy-analogy approximation', is equivalent 
to the exact solution of the Falicov-Kimball model (^0[). The Green function 
in the (full) Hubbard-III approximation includes the 'resonance broadening 
corrections' which lead to the following cubic equation for the Green function 

3[Gu-m{uj)f-8u[Gn-m{uj)?+'lGu-miu;){3+u;^-Uy4)-8Lu = . (105) 

This approximation becomes exact to (l/U)^, i.e., Gh-iii(^;) —>■ Go{z) for 
U — > oo, cf. (|2q). However, the Hubbard-III approximation predicts the 
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Fig. 3. Mott-Hubbard gap for one-particle excitations in the Hubbard 
model. The result from the expansion to second order in 1/C/ (solid line), 
from the Hubbard-III approximation (dotted line), and from the local- 
moment approximation (dashed line) are shown. 

transition to occur at C^"""^ = 2^/3 = 3.46 = 0.866VF. Although this is an 
improvement over the alloy- analogy approximation, the Hubbard bands still 
tend to repel each other, in contrast to our analytical result ( |lOl ) to 0{U-^). 

The one-particle gap in the Hubbard model as a function of U is shown 
in Fig. m for the Hubbard-III approximation and for the 1/U expansion. In 
the Hubbard-III approximation, the gap has a positive correction in Ai(f7), 

AH-in([/) = C/-4+i^, (106) 

in contrast to the exact first-order result (|8^). The Hubbard-III approxima- 
tion overestimates the stability of the Mott-Hubbard insulator already to 
first order in 1/U so that [4""^" < Uc appears to be natural. 

In Fig. ^ we compare the density of states for the Hubbard-III approx- 
imation and the 1/U expansion for U = 6 = {3/2)W. It is seen that the 
Hubbard-III approximation not only underestimates the gap, it also fails 
to reproduce qualitatively the tilt in the overall structure of the density of 
states. 
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1/U - Expansion 

Hubbard III - Approx. 




Fig. 4. Density of states in the Hubbard model for U = 6. The result 
from the expansion to second order (solid line) and from the Hubbard-III 
approximation (dashed line) are shown. 



6.2. Local-Moment Approach 

As a second example we compare the results of our systematic 1/U ex- 
pansion to those of the local-moment approach@, which provides an excel- 
lent description of the dynamical properties of the single-impurity Anderson 
model, particularly for strong couplingE3. 

The gap in LMA acquires the valuesHS 

^LMA^jj ^ 11 313) ^ 7 ^ A^^^{U = 8.4658) = 4.226 (107) 
so that we may approximate 

A^^HU) = U-A-'-f-'^. (108) 



The comparison with ( |103| ) shows that the first-order correction is very close 
to the exact first-order result. 

The one-particle gap as a function of U is shown in Fig. |3|. The neg- 
ative second-order contribution in ( |108D makes the local-moment gap bend 
downwards from the second-order result ( |103| ). Therefore, the local- moment 
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Fig. 5. Density of states in the Hubbard model for U = 4\/2 ~ 5.66. The 
result from the expansion to second order (dashed line) and from the local- 
moment approach (solid line) are shown. 



approach predicts the collapse of the insulator to occur at U^^^ 4.82 
1.21W which is larger than our second-order prediction Uc = 4.24 1.06H^ 
in (104). Since the Hubbard-III approximation underestimates Uc and the 
local-moment overestimates Uc, we are confident that the exact critical in- 
teraction strength cannot be much different from Uc ■ 

In Fig. |5| we display the density of states for the local-moment insulator^ 
for U = 4\/2 ~ 5.66 ~ lAW together with the corresponding result (|101| ) of 
the 1/U expansion to second order. Note that the energy scale is t* = \plt 
in Fig. m, and the density of states is rescaled accordinglyil. It is seen that 
the overall agreement is very good. The deviations are largest in the vicinity 
of the band edges, i.e., the single-particle gap deviates more than the overall 
density of states. In fact, the two curves are almost indistinguishable for U = 
6\/2 ~ 8.49 ~ 2.1W. This shows that the local-moment approach provides 
a very reliable description of the Mott-Hubbard insulator for U > lAW. 
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7. CONCLUSIONS 

In this work we have formulated the 1/U expansion for the single- 
particle Green function for the Mott-Hubbard insulator on the Bethe lattice 
with infinite coordination number. For this case the terms to order l/f/"" 
in the Kato-Takahashi perturbation theory can be replaced by two nth- 
order polynomials in /iq, one for the shape-correction correction terms and 
one for the gap-renormalization terms. The latter polynomial has only even 
(odd) powers for n odd (even). Therefore, the density of states, which is an 
important ingredient for the dynamical mean-field theory, is actually char- 
acterized by a (small) set of numbers for U ^ Uq- This observation should 
be useful for an assessment, and possible improvement, of the quality of nu- 
merical techmnueswhere the continuous density of states is replaced by a 
few peakHSSiS. 

We have tested our approach against the exact solution for the Falicov- 
Kimball model. The 1/U expansion up to second order gives an excellent 
description for U > 3t = 1.5U^^. The extrapolated gap is only 10% off the 
exact result. We are confident that the 1/U expansion is equally reliable 
for the Hubbard model in the regime Uq/U > 2/3, where we estimate Uc = 
(1.05 lb 0.05) from the extrapolation of the gap to second order in 1/U. 

As a first application of our result, we have tested two approximate the- 
ories for the Mott-Hubbard insulator which are exact to order (l/U)^. The 
first-order correction to the gap in the Hubbard-HI approximation resembles 
that of the Falicov-Kimball model, i.e., it has the wrong sign. Therefore, the 
density of states in the Hubbard-HI approximation qualitatively disagrees 
with the 1/U expansion, and the value for the Mott-Hubbard transition, 
C/H-in = 2^/3W ^ 0.866VF is too small. 

The density of states from the local-moment approach agrees very well 
with our results down to C/ ~ lAW . The first-order corrections to the gap 
are very close but the local-moment approach has a sizable second-order 
correction which is actually absent in the Hubbard model. Therefore, the 
local-moment estimate for the Mott-Hubbard transition is somewhat too 
large, U^^^ ^ 1.21W . 

There are two applications of our results. First, our formulae can be 
used as an input guess for the iteration cycle of the dynamical mean-field 
theory. This may help to stabilize and speed up the convergence of the it- 
eration in the Mott-Hubbard insulator for all U > Uc- Second, and more 
importantly, approximate theories on the Mott-Hubbard transition have to 
pass the test against our results for the single-particle density of states in the 
Mott-Hubbard insulator. A more detailed comparison of the existing nu- 
merical and analytical approximations is planned to be published elsewhere. 
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